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Abstract 

We study the time evolution for the quantum harmonic oscillator subjected to a sudden change 
of frequency. It is based on an analytic solution to the time dependent Ermakov equation for an 
approximated step function, and therefore, ours is a continuous treatment that differs from former 
studies that have limiting procedures for times when the step occurs. 
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I. INTRODUCTION 



The problem of the harmonic oscillator with time dependent frequency has received con- 
siderable attention over the years [1]. In particular Kiss et al. [2] have studied the problem 
of the time dependent frequency given by a step function (see also [3]). They have solved the 
problem by dividing it in two regions for the two different frequencies of the step function 
and matching later the solutions for the two regions. Here we would like to study the same 
problem but using a continuous approach, by using instead an invariant formalism. This will 
allow us to obtain analytic solutions that show squeezing (departing from an initial coherent 
state) or not depending on time. We will also show that, in a transformed Hilbert space, an 
initial coherent state remains coherent during the evolution. 

II. COHERENT STATES FOR CONSTANT FREQUENCY 

Let us first consider the Hamiltonian for the harmonic oscillator with unitary mass, 
m — 1, and unitary frequency, ujq — 1 

H^\{e+f). (1) 

We can define annihilation and creation operators as 

1 - 1 

h = -jJ,q + ip), = -j={q-ip), (2) 

such that we can rewrite the Hamiltonian as (we set h= 1) 

H^u;o(h^b+^^ =o;o(^n+0 . (3) 

Eigenstates for the Hamiltonian (3) are called Fock or number states: 

H\n) ^ LUo ( n + ^) \n) . (4) 



2^ 

Fock states are orthonormal and form a complete basis, such that any other state of the 
harmonic oscillator may be written in terms of them. In particular coherent states may be 
written as 

2 °° n 
n=0 * 
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where 

D{a) = e'*^^-"'"^ (6) 
Coherent states are eigenstates of the annihilation operator 

b\a) = a\a). (7) 

and have the property that the motion of the center of mass of the wave packet obeys the 
laws of classical mechanics (see for instance 



{a\q\a)=qc, {a\p\a) = Pc, {a\H\a) ^ (8) 
where the index c labels the classical variables. 

III. TIME DEPENDENT HARMONIC HAMILTONIAN 

The time dependent harmonic Hamiltonian reads 

H^l[f + n^{t)f]. (9) 
It is well-known that an invariant may be written for this interaction that has the form 



^ 1 
7= - 
2 



2 

<1\ , ...^2 



- ) +{pp- pqY 



(10) 



where p obeys the Ermakov equation 

p + n^p^p-\ (11) 

Furthermore, it is easy to show that / may be related to the Hamiltonian (1) by a unitary 
transformation of the form 

f = (12) 

that can be re-written as (see appendix) 
as 

E = (jJofifK (14) 



By using Eq. (4) we can see that 



if^n)^u;o(n+^^ ft 



n) 



(15) 



i.e. states of the form 

\n)t^f^n) (16) 

are eigenstates of the so-called Lewis invariant. Lewis [4] wrote his invariant in terms of 
annihilation and creation operators 



/ = a/ a + - = + t:- 



(17) 



with 



V2 



^ + i{pp- pq) 



V2 



-^-i{pp- pq) 



(18) 



Once defined creation and annihilation operators, analogous equations to the harmonic 
oscillator with constant frequency may be obtained for the TDHO, for instance 



Dtia) 



aa^ —a*a 



(19) 



and 



with 



a\a)t — Oi\a)t- 



\a)t = e 2 



n \ Tl. 



n=0 



(20) 



(21) 



Recently it has been shown that the Schrodinger equation for the time dependent har- 
monic Hamiltonian has a solution of the form [5] 



|^(i)) = e-^^"-''o'^W'^*f'tj'(o)|^(0)). 
with ijj{t) — 1/ p^. If we consider the initial state to be 

|V'(0))=ft(0)|a) = |«)o 
with a given in (5), we note that the evolved state has the form 



(22) 



(23) 



(24) 



this is, coherent states keep their form through evolution. 
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IV. MINIMUM UNCERTAINTY STATES 



We define the operators 

g = i=(a + at), (25) 

and 

P=^(a-a.). (26) 
It is easy to see that they are related to q and p by the transformations 

q = fgr^ (27) 

and 

p = fPf^ (28) 

and that they obey the equations [7] 

P = iuj(t)[i, P] = -uj{t)Q, (29) 

and 

Q^iuj{t)[i,Q]=uj{t)P. (30) 
The uncertainty relation for operators Q and P for the coherent state (21) is given by 

AQAP = i (31) 



where AX = Y (X^) — (Time dependent) Coherent states are thus minimum uncer- 

tainty states (MUS), not for position and momentum but for the transformed position and 
momentum. 



V. STEP FUNCTION 



Let us consider the Hamiltonian of the system to be given by equation (9) with Q{t) 
given by a step function that may be modeled by [6] 



n{t) = uji 



l + ^(l + ^tanh[e{t-ts)] 



(32) 



where tg is the time at which the frequency is changed, A — 002 — ^^i with cui and CU2 the 
initial and final frequencies, and e is a parameter (e — > 00 models better the step function). 



In fig. (1) we plot this function as a function of t for cui — 1 and a;2 = 2 (solid line) and 
u;2 — S (dash line). The solution to the Ermakov equation (11) for this particular form of 
Q{t) is given by [6] 




A plot of p{t) is given in fig. 2 for the same values as Fig. (1). We also plot uj{t) = l/p^ in 
fig. (3). It may be numerically shown that the time average of omega{t) from ts to the end 
of the first period is 2 for the solid line and 3 for the dash line, with cumax — 2^ for the solid 
line and cUmax — 3^ for the dash line. 

Let us consider that at time t — we have the system in the initial state (5). Prom 2 
we can see that f{0) = 1 as p = and Inp = 0. Therefore from (23) |V'(0)) = ^^(O)!^) = 
|q;)o = I a) and from (22) we obtain the evolved wave function 

|^(t)) = e-'^^o^^t)dtf^\a) = rt|cte-^^o'^W<i*). (34) 

Note that the coherent state in the above equation is given in the original Hilbert space, i.e. 
in terms of number states given in (4). Prom fig. 2 we can also see that for the maximums, 

p{tmax) = and liap{tmax) = 0, therefore f^{tmax) = 1 and 

mt^ax)) = \ae-^^o-^m^^ (35) 

i.e. we obtain back the initial coherent state. However, for the minimums, we have p{tmin) — 
and \n.p{tmin) 7^ and then we obtain 

mmin)) = e'^^'^(^^+^^)|ae--^o-^"-W'^*), (36) 

that may be written in terms of annihilation and creation operators as 

mmin)) = e^^'^(^^-(^')')|ae-/o-™-W<^*), (37) 

that are the well-known squeezed (two-photon coherent) states [10] 

|^(t.n.n)) = |ae-/o-^"-W'^*;lnp(t^,J) (38) 

that, as coherent states, squeezed states are also MUS. However the uncertainties for q and 
p are 
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V2p{tmin) 

and 

Ap = (40) 
v2 

For times in between we will have neither coherent states nor standard squeezed states (in 
the initial Hilbert space), but the wave function 

It should be stressed however, that in the instantaneous Hilbert space we will always have 
the coherent state (24). 



VI. CONCLUSIONS 

We have studied the problem of the time dependent harmonic oscillator for a particular 
form of time dependency, namely the step function. We have studied it from an invariant 
point of view that has made it possible to obtain analytic solutions. We have shown that, 
depending on the space we look at the solutions for an initial coherent state, the state 
remains coherent or it may present squeezing. This squeezing may be enhanced by making 
larger the change in frequency. 

VII. APPENDIX 

The unitary operator T is given by 

f = p^^i^P+P^+^J') ^ /-¥i^P+P^+^J')_ (42) 

The exponential above has the sum of the operators 

A^i^^{qp + pq) (43) 



and 

B = tlnp-^^f. (44) 

One can show that 

[A,B]^2\npB. (45) 



Using the fact that the commutator of the two operators is proportional to one of the 
operators, we can factor the exponential in the form (see for instance [9]) 



FIG. 1: ri(t) as a function of t for = 1 and a;2 = 2 (solid line) and 102 = 2i (dash line), e = 20 
and tg = 2. 

FIG. 2: p{t) as a function of t for uji = 1 and u)2 = 2 (solid line) and ^2 = 3 (dash line), e = 20 
and tg = 2. 

FIG. 3: w(t) as a function of t for wi = 1 and = 2 (solid line) and W2 = 3 (dash line), e = 20 
and tg = 2. 
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